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1 Introduction 



Heavy- ion collisions (HIC) provide a unique opportunity to study the properties of hot 
and dense nuclear matter as the nuclear matter is produced at temperature, densities 
and neutron-proton ratios away from that of the ground state nuclei. The nuclear equa- 
tion of state (EOS) as well as the in-medium nucleon-nucleon (nn) cross section plays a 
great role in determining these properties [1,2]. The importance of collective transverse 
in-plane flow in determining these properties has been known for past three decades. 
The collective transverse in-plane flow has been found to be influenced by the entrance 
channel parameters like mass of the colliding system [3,4], impact parameter [5-7] and 
incident energy [7-9] . The competition between attractive interactions due to mean fleld 
(at low incident energies) and repulsive interactions due to nucleon-nucleon (nn) colli- 
sions (at high incident energies) led to the disappearance of flow at a particular incident 
energy (where attractive and repulsive interactions counterbalance each other), known 
as energy of vanishing flow (EVP) [10]. Since one of the motivations behind studying 
HIC is to extract information about the hot and dense nuclear matter, so few studies 
also shed light on the properties such as density and temperature reached in HIC and 
also on the participant-spectator matter [11, 12]. Puri and collaborators have shown that 
participant-spectator matter can act as a barometer for the study of EVP. The study 
also pointed towards the insensitivity of participant-spectator matter towards equation 
of state, momentum-dependent interactions (MDI) and nn cross sections. These studies 
were done at the EVP. Moreover, these studies were restricted to the systems lying close 
to the stability line. 

With the existing and upcoming radioactive ion beam (RIB) facilities around the 
world, the present physics interest have been shifted towards nuclei far from the stability 
line. Various phenomena like fusion-flssion, collective flow, EVP and multifragmenta- 
tion have been studied for such nuclei which lie far from the stability line. So, in the 
present paper, we aim to see the properties like, density, temperature, thermalization, 
and participant-spectator matter for neutron-rich systems lying far from the stability 
line. 

The present study is carried out within the framework of isospin-dependent quantum 
molecular dynamics (IQMD) model which is described in detail in section 2. Our re- 
sults along with discussions are presented in section 3 and we summarized the results in 
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section4. 



2 The model 

The IQMD model [13] which is the extension of quantum molecular dynamics (QMD) 
[14] model treats different charge states of nucleons, deltas, and pions explicitly, as in- 
herited from the Vlasov-Uehling-Uhlenbeck (VUU) model. The IQMD model has been 
used successfully for the analysis of a large number of observables from low to relativistic 
energies. Puri and coworkers have demonstrated that QMD and IQMD carries essential 
physics needed to demonstrate the various phenomena such as collective flow, multifrag- 
mentation and particle production [15, 16]. The isospin degree of freedom enters into the 
calculations via symmetry potential, cross sections and Coulomb interactions. 

In this model, baryons are represented by Gaussian-shaped density distributions 

1 1 2L 

fi{r,P,t) = l^^^'M-lr-nit)]^—) X exp(-[p-pl(t)]2— ) (1) 

where L is the Gaussian width which can be regarded as a description of the interaction 
range of a particle. In IQMD the Gaussian width is system dependent. The system 
dependence of L in IQMD has been introduced in order to obtain the maximum stability 
of the nucleonic density distribution, for Au-|-Au a value of L = 8.66 fm2 is chosen, for 
Ca+Ca and lighter nuclei L = 4.33 fm2. Nucleons are initialized in a sphere with radius 
R = 1.12 A^/^ fm, in accordance with liquid-drop model. Each nucleon occupies a volume 
of h^, so that phase space is uniformly filled. The initial momenta are randomly chosen 
between and Fermi momentum (pV). The nucleons of the target and projectile interact 
by two- and three-body Skyrme forces, Yukawa potential and Coulomb interactions and 
momentum-dependent interactions. In addition to the use of explicit charge states of all 
baryons and mesons, a symmetry potential between protons and neutrons corresponding 
to the Bethe-Weizsacker mass formula has been included. The hadrons propagate using 
Hamilton equations of motion: 

dfj _ d{H) _ dpi _ d{H) 
dt dpi ' dt dfi 

with 
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(H) = {T) + {V) 

i i j>i 

xfj{f',p',t)dfdf'dpdp'. (3) 



The baryon potential V'-" in the above relation, reads as 



Y^3 r) — y Skyrme + ^Yukawa + ^Coul + ^mdi + Kyi 

= [ti5{r'-r)+t26{r'-r)p''-\ — 
exp(|(f'-f)|//i) , Z,Z,e^ 



sym 



(l(f'-f)l^) ' |(f'-f)| 

+U\Ti^[h{p' -p)^ + l]5{f' -r) 

+te-TsinjS{fi' -rj). (4) 

The Skyrme potential y^skyrme consists of two-body and three-body interactions and can 
be parameterized as: 



U = a 



PoJ \PoJ 



The parameters a and ^ can be fixed by requirement that at normal nuclear matter 
density the average binding energy should be -16 McV and total energy should have a 
minimum at po- Different values of 7 lead to different equation of states. Naturally, 
larger value of 7 leads to hard EOS whereas smaller value of 7 results in soft equation of 
state. The compressibility K (MeV) is 200 and 380 for soft and hard EOS, respectively. 
Here te = 100 MeV and Zi and Zj denote the charges of ith and jth baryon, and T^i 
and Tsj are their respective components (i.e., 1/2 for protons and —1/2 for neutrons). 
The parameters /i and ti,....,tQ are adjusted to the real part of the nucleonic optical 
potential. The Yukawa potential Vyukawa IQMD is very short ranged (/x = 0.4 fm) 
and weak. For the density dependence of the nucleon optical potential, standard Skyrme 
type parametrization is employed. We use a soft equation of state with and without 
momentum-dependent interactions (MDl) labeled as Soft and SMD, respectively along 
with the standard isospin- and energy-dependent nn cross section reduced by 20%, i.e. a ~ 
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0.8 (T^^^. The cross section for proton-neutron collisions is three times as that for proton- 
proton or neutron-neutron collisions. The cross sections for neutron-neutron collisions 
are assumed to be equal to the proton-proton cross sections. The more details about the 
elastic and inelastic cross sections for proton-proton and proton-neutron collisions can be 
found in [13,17]. In a recent study, Gautam et al. [18] has confronted the theoretical 
calculations of IQMD with the data of ^^Ni Ni and ^^Fe Fe [19]. The results 
with the soft EOS (along with the momentum-dependent interactions) and above choice 
of cross section are in good agreement with the data at all colliding geometries. Two 
particles collide if their minimum distance d fulfills 



where 'type' denotes the ingoing coUision partners (N-N....). Exphcit Pauh blocking is 
also included; i.e. Pauh blocking of the neutrons and protons is treated separately. We 
assume that each nucleon occupies a sphere in coordinate and momentum space. This 
trick yields the same Pauli blocking ratio as an exact calculation of the overlap of the 
Gaussians will yield. We calculate the fractions Pi and P2 of final phase space for each 
of the two scattering partners that are already occupied by other nucleons with the same 
isospin as that of scattered ones. The collision is blocked with the probability 



and, correspondingly is allowed with the probability 1 - ^ block- For a nucleus in its 
ground state, we obtain an averaged blocking probability {Puock) = 0.96. Whenever an 
attempted collision is blocked, the scattering partners maintain the original momenta 
prior to scattering. 

3 Results and discussion 

We simulate the reactions of Ca+Ca and Xe-|-Xe series having N/Z ratio varying from 1.0 
to 2.0 in steps of 0.2. In particular we simulate the reactions of ^^Ca-l-^^Ca, "^**Ca-|-"^*^Ca, 
48Ca+48Ca, 52Ca+52Ca, seCa+^eCa and 

60Ca+60Ca and i^oXe+^^oXe, ^^oxe+i^oxe, i29xe+i29xe, i^oxe+i^oxe, i^iXe+i^iXe, and 
^^^Xe-F^^^Xe at h/\>max = 0.2-0.4. The incident energies are taken to be 50 and 250 
Me V/ nucleon. 




(6) 



Phiock = 1 - [1 - min{Pi, 1)][1 - min{P2, 1)], 



(7) 
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Figure 1: (Color online) The N/Z dependence of participant and spectator matter for the 
reactions of Ca+Ca at 50 (upper panels) and 250 (lower) MeV/nucleon for different N/Z 
ratios with Egym = MeV (left panels) and MDI (right). Various symbols are explained 
in the text. 



In the present work, we define the participant and spectator matter in terms of the 
nucleonic concept. All the nucleons which have experienced at least one collision are la- 
beled as participant matter, scaled to the total mass of the reacting nuclei. The remaining 
matter is labeled as spectator matter. These definitions provide a possibility to study the 
reaction in terms of participant-spectator fireball model. 

In fig. 1, we display the N/Z dependence of participant-spectator matter for Ca+Ca 
reactions. Upper (lower) panels represent the results for 50 (250) MeV/nucleon. Solid 
(open) symbols are for participant (spectator) matter. From figure (circles), we find that 
at lower incident energy there is a slight change in participant-spectator matter with N/Z 
ratio whereas, at higher incident energy of 250 MeV/ 

nucleon, the participant-spectator matter is almost independent of the neutron content 
of the colliding systems. As we are increasing the neutron content at fixed Z, the mass of 
the system also increases. Due to increase in mass, the participant matter should increase 
with N/Z ratio (at fixed incident energy). But as we are moving to higher neutron 
content, the role of symmetry energy also increases, which being repulsive in nature leads 
to lesser number of collisions in systems with more neutrons and hence participant matter 
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should decrease with increase in neutron content of the system. So the net effect is due 
to the interplay of mass effects and effect of the symmetry energy. To check; the relative 
importance of these two effects we make the strength of potential part of symmetry energy 
zero and calculate the participant-spectator matter for Ca+Ca reactions. The results are 
displayed in Fig. 1(a) (50 MeV/nucleon) and 1(b) (250 MeV/nucleon) (triangles). Prom 
figure, we see that on reducing the strength of symmetry energy to zero, the participant 
(spectator) matter increases (decreases) slightly for higher neutron content whereas the 
effect on lower N/Z ratios is almost negligible. 

The momentum-dependent interactions also play a significant role in the dynamics of 
heavy-ion collisions. To see the role of MDl on participant-spectator matter we calculate 
the participant-spectator matter for Soft-|-MDI (SMD) EOS for the reactions of Ca-|-Ca. 
The results are displayed by diamonds in Fig. 1(c) and 1(d) at 50 and 250 MeV/nucleon, 
respectively. From figure, we see that participant [closed diamonds] (spectator, open 
diamonds) matter decreases (increases) with SMD. This is due to the fact that since MDI 
is repulsive in nature, so the matter is thrown away from the central dense zone and hence 
lesser number of collisions will take place which leads to decrease (increase) in participant 
(spectator) matter. Similar behaviour is observed for Xe-|-Xe reactions (results not shown 
here) . 

Since one of the motivations behind studying heavy-ion collisions is to extract in- 
formation about the hot and dense nuclear matter. In our approach, matter density is 
calculated by [20] 

''(f.*) = T'p^e(-('--'-.('))V-). (8) 

Here At and Ap stand, respectively, for the mass of the target and projectile. In actual 
calculations, we take a sphere of 2 fm radius around the center-of-mass and compute the 
density at each time step during the reaction using Eq. (5.1). Naturally, one can either 
extract an average density (p""^) over the whole sphere or a maximal value of the density 
l^pTuax^ reached anywhere in the sphere. 

In fig. 2, we display the time evolution of maximum (upper panels) and average 
density (lower panels) at 50 (left panels) and 250 MeV/nucleon (right) for the reactionsof 
Ca-|-Ca. Solid, dashed, dotted, dash-dotted, short dotted and short dash-dotted represent 
N/Z ratios of 1.0, 1.2, 1.4, 1.6, 1.8, and 2.0, respectively. From fig. 2 (left panels), we 
find that maximum and average density rises as the reaction proceeds, reaches maximum 
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Figure 2: (Color online) The time evolution of the maximum (< Pmax/ Po >) and average 
(< pmavgl >) densities for the reactions of Ca+Ca at 50 (left panels) and 250 (right) 
MeV/nucleon. Lines have been explained in text. 

at 10-20 fm/c when the matter is highly compressed and finally decreases during the 
expansion phase. The average density reached in the reactions is about 1.5 times the 
normal nuclear matter density. The maximal value of the density decreases with increase 
in neutron content of the system. This slight decrease of the maximal value may be due 
to the larger effect of symmetry energy in systems with higher neutron content, which 
in turn, pushes the matter away from the central dense zone due to the repulsive nature 
of symmetry energy. From fig. 2 (right panels), we see that at 250 MeV/nucleon, much 
higher density is achieved. This is because of the much violent nature of the reaction at 
higher incident energy. 

In fig. 3, we display the N/Z dependence of the maximal value of the maximum and 
average density achieved in a reaction at 50 (upper panels) and 250 (lower) MeV/nucleon 
for the reactions of Ca+Ca reactions. Solid (open) circles represent the maximal value 
of Pavg (Pmax)- Wc find that maximal value almost remains independent of the neutron 
content of the reacting system. 

To see the role of MDI and symmetry energy, we display the results with MDI (Figs. 
3(a) and 3(b)) (diamonds) and with no symmetry energy (Figs. 3(c) and 3(d)) (triangles) 
on average (closed symbols) and maximum density (open symbols) achieved in reactions 
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Figure 3: (Color online) The N/Z dependence of maximal value of maximum (< 
Pmax/po >) and average density (< Pavg/po >) for the reactions of Ca+Ca at 50 (up- 
per panels) and 250 (lower) MeV/nucleon with MDI (left panels) and Egym = MeV 
(right). Symbols have same meaning as in Fig. 1. 
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Figure 4: (Color online) The N/Z dependence of maximal value of maximum (Tm,aa;)and 
average temperature (Tavg) for Ca+Ca (upper panels) and Xe+Xe (lower) at 50 (left) 
and 250 (right) MeV/nucleon. Diamonds and triangles represent the results for MDI and 
Esj/m = MeV, respectively. 
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of Ca+ 

Ca. From figure, we see that average and maximum density decreases with SMD as 
compared to soft EOS, although the N/Z behavior remains independent. Similarly, when 
we reduce the strength of symmetry energy to zero, the average and maximum density 
increases slightly for higher neutron content whereas for lower neutron content, the effect 
of symmetry energy is negligible. 

The associated quantity linked with the dense matter is the temperature. In principle, 
a true temperature can be defined only for a thermalized and equilibrated matter. Since 
in heavy-ion collisions the matter is non-equilibrated, one can not talk of "temperature" . 
One can, however, look in terms of the local environment only. In our present case, we 
follow the description of the temperature given in Refs. [20,21]. In the present case, 
extraction of the temperature T is based on the local density approximation, i.e., one 
deduces the temperature in a volume element surrounding the position of each particle at 
a given time step [20, 21]. Here, we postulate that each local volume element of nuclear 
matter in coordinate space and time has some "temperature" defined by the diffused edge 
of the deformed Fermi distribution consisting of two colliding Fermi spheres, which is 
typical for a nonequilibrium momentum distribution in heavy-ion collisions. 

In this formalism (dubbed the hot Thomas- Fermi approach [20]), one determines ex- 
tensive quantities like the density and kinetic energy as well as entropy with the help of 
momentum distributions at a given temperature. Using this formalism, we also extracted 
the average and maximum temperature within a central sphere of 2 fm radius as described 
in the case of density. 

In fig. 4, we display the N/Z dependence of the temperature achieved in a reaction. Upper 
(lower) panels display the results for Ca-|-Ca (Xe-|-Xe) reactions. Solid (open) symbols 
represent the results of T""^ (^maxy j^g^^ (right) panels represent the results of 50 (250) 
MeV/nucleon. We find that higher temperatures are achieved in reactions taking place at 
higher incident energy. For low incident energy, temperature remains almost constant for 
lower neutron content and shows slight increase for higher neutron content. We also study 
the effect of MDI (diamonds) and symmetry energy (triangles) on average temperature 
for the reactions of Ca+Ca at 50 and 250 MeV/nucleon. From figure, we see that aver- 
age temperature decreases for all N/Z ratios when MDI is included. Similarly, when we 
reduce the strength of symmetry energy to zero, there is a very slight decrease in average 
temperature. The N/Z behavior, although, remains unaltered. Similar behavior of MDI 
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Figure 5: (Color online) The time evolution of anisotropy ratio for Ca+Ca (upper panels) 
and Xe+Xe (lower) reactions at 50 (left panels) and 250 (right) MeV/nucleon for different 
N/Z ratios. Lines have same meaning as in Fig. 2. 



and symmetry energy is observed for Xe+Xe reactions (results not shown here). 

In fig. 5 we display the time evolution of anisotropy ratio < Ra > reached in a reaction. 
The < i?a > is defined as 



2\ . . 

This anisotropy ratio is an indicator of the global equilibrium of the system. This 
represents the equilibrium of the whole system and does not depend on the local positions. 
The full global equilibrium averaged over large number of events will correspond to < Ra > 
= 1. From figure, we see that as the reaction proceeds, the < Ra > ratio increases and 
saturates after the high dense phase is over. We also find that the effect of N/Z ratio is 
very less on the thermalization achieved in a reaction. We also find a smaller value of 

< -Ra > at the start of the reaction for higher energy projectile (see right panels). 

In fig. 6, we display the time evolution of relative momentum. The relative momentum 

< Kn > of two colliding Fermi spheres, is defined as 

< Kr >=< \Pp{f,t) - PT{r,t)\/n >, (10) 
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Figure 6: (Color online) Same as Fig. 5, but for the time evolution of relative momentum. 



where 

P,{f,t) = ^'-' ^ = l,2, (11) 

Pj{r,t) 

i.e., i = 1 and 2 correspond to ^ = P (Projectile) and T(target). Here Pj and pj are the 
momentum and density of the jth particle and i stands for either projectile or target. 
The <Kfj> is an indicator of the local equilibrium because it depends also on the local 
position r. 

We find that relative momentum decreases as the reaction proceeds. The effect of 
neutron content is very less on the value of relative momentum achieved in a reaction. 
The effect diminishes further at higher incident energy. Also, at higher incident energy, 
at the start of the reaction we have a higher value of the relative momentum. Also the 
decrease in <K^> is much faster for higher incident energy. 

4 Summary 

We studied the dependence of participant-spectator matter, temperature and density 
reached in heavy-ion reactions of neutron-rich systems on the neutron content of the 
colliding system. The study is done at 50 and 250 MeV/nculeon. Our calculations showed 
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that the behaviour of these quantities is almost independent of the neutron content of the 
system under study. In addition, we also studied the thermalization achieved in neutron- 
rich systems. Thermalization also showed insensitivity towards the neutron content of 
the reacting systems. 

This work has been supported by a grant from Centre of Scientific and Industrial 
Research (CSIR), Govt, of India. Author is thankful to Profs. J. Aichehn and R. K. Puri 
for enlightening discussions on the present work. 
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